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Abstract: Motivated by our earlier paper [Q, we discuss how the infamous grav- 
itino problem has a natural built in solution within supersymmetry. Supersymmetry 
allows a large number of flat directions made up of gauge invariant combinations of 
squarks and sleptons. Out of many at least one generically obtains a large vacuum 
expectation value during inflation. Gauge bosons and Gauginos then obtain large 
masses by virtue of the Higgs mechanism. This makes the rate of thermalization 
after the end of inflation very small and as a result the Universe enters a quasi- 
thermal phase after the inflaton has completely decayed. A full thermal equilibrium 
is generically established much later on when the flat direction expectation value has 
substantially decareased. This results in low reheat temperatures, i.e., Tr ~ O(TeV), 
which are compatible with the stringent bounds arising from the big bang nucleosyn¬ 
thesis. There are two very important implications: the production of gravitinos and 
generation of a baryonic asymmetry via leptogenesis during the quasi-thermal phase. 
In both the cases the abundances depend not only on an effective temperature of the 
quasi-thermal phase ( which could be higher, i.e., T Tr), but also on the state 
of equilibrium in the reheat plasma. We show that there is no “thermal gravitino 
problem” at all within supersymmetry and we stress on a need of a new paradigm 
based on a “quasi-thermal leptogenesis”, because in the bulk of the parameter space 
the old thermal leptogenesis cannot account for the observed baryon asymmetry. 




Contents 


[l]. Introduction 


Thermalization in supersymmetric theories: Part-1 

A brief recourse to thermalization 


2 . 


2.2 Flat directions in Supersymmetry 


2.3 Flat directions during and after inflation 


2~4| Flat directions and inflaton decay 


2.5 Flat directions and thermalization 


0 

0 


n 

T2 

13 


B 

i 


Reheat temperature of the Universe 


3.1 Numerical examples 


Thermalization in supersymmetric theories: Part-2 

H Higher order superpotential terms and Kahler corrections 
|4.2| Models with gauge-mediated supersymmetry breaking 
PI Split supersymmetry 


|5j. Particle production during thermalization 
|5.1| Quasi-adiabatic evolution of the Universe 
|0 Quasi-thermal particle production 


T6 

17 


18 

T9 

2C 

22 


24 

24 

1261 


||. Gravitino production in a quasi-thermal case 

[?. Leptogenesis in a quasi-thermal case 

0 Basic concept 
|0 Quasi-thermal leptogenesis 


26 


m 

n 

|33| 


|8j. Thermalization after preheating 


35 


|9. Conclusion 


36 


10. Acknowledgments 


38 


1 






































































Tlj. Appendix 

m Last stage of inflaton decay 
11. 2| Gravitationally decaying inflaton 


11. 3| Right-handed sneutrino as the inflaton 


11.4j Flat directions and thermalization: additional considerations 


39 

39 

42 

1 

44 


1. Introduction 

Primordial inflation |2| 1 is the best paradigm for explaining the initial conditions for 
the structure formation and the cosmic microwave background (CMB) anisotropies Q. 
The Universe is cold and empty after inflation and all the energy is stored in the 
inflaton condensate oscillating around the minimum of its potential. Reheating is 
an important stage for any inflationary model. It describes the transition from this 
frozen state to a hot thermal Universe. It involves various processes and eventually 
results in a thermal bath of elementary particles which contains the Standard Model 
(SM) degrees of freedom. 

The Inflaton decay is the most relevant part of the reheating. Only one-particle 
decay of the non-relativistic inflaton quanta were considered initially || 2 . The 
treatment is valid if the energy transfer to the fields which are coupled to the inflaton 
takes place over many inflaton oscillations. This requires that the inflaton couplings 
to the SM fields are sufficiently small. Usually it is assumed that the plasma reaches 
complete kinetic and chemical equilibrium immediately after all the inflaton quanta 
have decayed. The reheat temperature of the Universe, T R , is then quoted as: [|oJ 

T R ~0.5g^WM P T d , ( 1 . 1 ) 

where g* is the total number of relativistic degrees of freedom, T c j is the inflaton 
decay rate and Mp = 2.4 x 10 18 GeV is the reduced Planck mass. 

It was also realized that the coherent oscillations of the inflaton can create par¬ 
ticles non-perturbatively 01- This mechanism is called preheating and it is par¬ 
ticularly efficient when the final products are bosonic degrees of freedom. It only 

1 For a realistic toy model of inflation where the inflaton carries the charges of minimal super- 
symmetric Standard Model, and gives rise to a low scale inflation, see [||. 

2 For a recent discussion on supersymmetric reheating, see ||. 


2 





































takes about two dozens of oscillations to transfer the energy from the homogeneous 
condensate to non-zero modes of the final state(s) ||. However, despite efficiently 
transferring the energy, preheating does not result in a complete decay of the inflaton. 
An epoch of perturbative reheating is an essential ingredient of any potentially realis¬ 
tic cosmological model || 3 4 . In supersymmetry (SUSY) a first stage of preheating 
is naturally followed by the last stage of perturbative decay (see Appendix |1 1. 1|) . 
For these reasons we concentrate on the perturbative inflaton decay, because of its 
relevance to create a thermal bath of SM particles 3 4 5 . 


Often cosmology is considered as a probe to the early Universe, a well known ex¬ 
ample is the Gravitino Problem in the context of a supersymmetric cosmology. SUSY 
introduces new degrees of freedom and new parameters. Most of them are rather 
poorly constrained from experiments. Cosmology however acts as a test bed where 
some of the SUSY particles can be tested. In this regard the reheat temperature 
plays an important role as we shall explain below. 


The gravitino is a spin-3/2 supersymmetric partner of the graviton, which is 
coupled to the SM particles with a gravitational strength. The Gravitinos with both 
the hclicities can be produced from a thermal bath. There are many scattering 
channels which include fermion, sfermion, gauge and gaugino quanta all of which 
have a cross-section oc 1/M'f |14], which results in a gravitino abundance (up to a 


3 Preheating ends due to backreaction as well as the expansion of the universe. Preheating does 
not destroy the zero mode of the inflaton condensate completely, although the amplitude of the 
inflaton oscillations diminish, but the inflaton decay is completed when the zero mode perturbatively 
decays into the SM or some other degrees of freedom, see Under only special conditions the 
inflaton condensate can completely fragment, see |jTo| . 

4 Besides, the initial stages of preheating can give rise to a large non-Gaussianity which already 
rule out the bulk of the parameter space from the current observations jll), see for a detailed 
discussion in our recent paper Jl2| . 

5 In supergravity inflationary models the inflaton can be a hidden sector which interacts gravita¬ 
tionally with other sectors. In which case the inflaton decay is perturbative at all times. In string 
theory however the epoch of thermalization could be very complicated as described in the recent 
paper jli|. This is due to the fact that string thermodynamics differs in many respects to that of 
a particle thermodynamics. 
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where Mg is the gluino mass. Note that for Mg < m 3/ / 2 both the helicity states 
have essentially the same abundance, while for Mg m 3/ / 2 production of helicity 
±1/2 states is enhanced due to their Goldstino nature. The linear dependence of the 
gravitino abundance on T R can be understood qualitatively. Since the cross-section 
for the gravitino production is oc Mp 2 , the production rate at a temperature, T, 
and the abundance of the gravitinos produced within one Hubble time will be oc T 3 
and oc T respectively. This implies that the gravitino production is efficient at the 
highest temperature of the radiation-dominated phase of the Universe, i.e., T R . 

An unstable gravitino decays to particle-sparticle pairs, and its decay rate is given 
by r 3/2 — mgy 2 /4Mp Hl4j| . If m 3 / 2 < 50 TeV, the gravitinos decay during or after 
big bang nucleosynthesis (BBN) ]T7] , which can ruin its successful predictions for 
the primordial abundance of light elements |L8|. If the gravitinos decay radiatively, 
the most stringent bound, (n 3 / 2 /s) < 1CV 14 — 10~ 12 , arises for m 3 / 2 — 100 GeV — 1 


TeV 19 


On the other hand, much stronger bounds are derived if the gravitinos mainly 
decay through the hadronic modes. In particular, for a hadronic branching ratio ~ 1, 
and in the same mass range, (n 3 / 2 /s) < 10” 16 — 1CV 15 will be required [|2(] . 


For a radiatively decaying gravitino the tightest bound (n 3 / 2 /s) < 10 -14 arises 
when m 3 /2 — 100 GeV [19]. Following Eq. (|1.2|) the bound on reheat temperature 
becomes: 7 r < 10 10 GeV. This turns out to be a very stringent limit if the inflaton 
decay products immediately thermalize. In fact, for the inflaton mass m$ = 10 13 
GeV, it is at best marginally satisfied even for a gravitationally decaying inflaton 
(see an example given in an Appendix |1 1. 2|) . For a TeV gravitino which mainly 
decays into gluon-gluino pairs (allowed when m 3 / 2 > Mg) a much tighter bound 
(^3/2/s) < io~ 16 is obtained [ 2 (]] r which requires quite a low reheat temperature: 
T r < 10 6 GeV. 

The gravitino will be stable if it is the lightest supersymmetric particle (LSP), 

where A-parity is conserved. The gravitino abundance will in this case be constrained 

6 From now on we take <?* = 228.75 as in the case of minimal supersymmetric SM (MSSM). 
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by the dark matter limit, [||], f2 3 / 2 /7 2 < 0.129, leading to 


'^ 3/2 


< 5 x 10“ 10 


1 GeV 

m 3/2 


(1.3) 


For 777 . 3/2 < Mg, the helicity ±1/2 states dominate the total gravitino abundance. 
As an example, consider the case with a light gravitino, 7773/2 = 100 KeV, which can 
arise very naturally in gauge-mediated models |2l|. If Mg ~ 500 GeV, see Eq. (|1.2|) , a 
very severe constraint, Tr < 10 4 GeV, will be obtained on the reheat temperature '. 

As we have argued, gravitinos indirectly put constraints on the thermal history 
of the Universe, which is often known as the gravitino problem , in models with weak 
scale SUSY. The reheat temperature, Tr, plays a very important role here as it 
sets the thermal (and possibly the largest) contribution to the gravitino abundance. 
However besides the gravitino production, it has other implications for cosmology, for 
instance, thermal leptogenesis, thermal production of weakly interacting dark matter 
particles, etc. They are all directly or indirectly connected to the reheat temperature 
of the Universe. 

The most important worry is the central assumption that, the SM and or 
MSSM degrees of freedom thermalize instantly. In our opinion this is the key 
assumption which has not been questioned well enough in the literature. It is only 
very recently that various issues of thermalization have been considered in a NON- 
SUSY case carefully We intend to cover this lapse and describe thermalization 
process within SUSY. The subject stands on its own right because SUSY introduces 
new ingredients into the game, i.e., flat directions. We will elaborate on the role of 
flat directions in thermalization, which was first pointed out in our earlier paper [;!,]. 

SUSY along with gauge symmetry introduces gauge invariant flat directions. I 11 
any supersymmetric extension of the SM there are flat directions primarily made 
up of squarks, sleptons (SUSY partners of SM quarks and leptons) and Higgses. A 
flat direction in a cosmological context can take a large vacuum expectation value 
(VEV) by virtue of a shift symmetry. During inflation quantum fluctuations of any 
light field accumulate in a coherent state and its VEV makes a random walk with its 
variance growing linearly in time. The condensate becomes homogeneous on scales 
larger than the size of the Hubble radius with a growing VEV. However in many 
cases (e.g., minimal supergravity) shift symmetry is not protected which truncates 

7 Similar results were obtained in connection to the brane worlds, see for instance |2^], and for 
the discussions on reheating, see pq] and p4) . 
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the VEV to be as large as the four dimensional Planck scale, Mp. Moreover SUSY 
is broken and the flat directions obtain soft mass terms and A-terms. It is also 
possible that non-renormalizable terms arise in the superpotential after integrating 
out heavy degrees of freedom which are associated to a new physics at high scales. 
These contributions lift the flatness but still allow the VEV to be significantly large 


in a wide class of models, for a review see [26 


After inflation the VEVs of flat directions do not settle at the origin immediately. 
The cosmologically sliding VEV plays a crucial role in slowing down thermalization 
of the inflaton decay products as mentioned in Ref. [[j]]. The idea is very simple, if 
a flat direction develops a VEV, the SM gauge fields become massive, similar to the 
Higgs mechanism. The masses for the gauge bosons (and gauginos) are proportional 
to the VEV of the flat direction. This obviously breaks the charge and the color in 
the early Universe, but the VEV of the flat direction finally vanishes and therefore 
there is no threat to the low energy phenomenology. 

Consider a situation where the inflaton has completely decayed into the MSSM 
degrees of freedom 8 . The resulting plasma is initially far away from a full thermal 
equilibrium. The leading reactions which establish equilibrium are the 2 —> 2 and 

The former lead to kinetic 


2 —> 3 scatterings mediated by gauge fields |25], |2£ 
equilibrium, while the latter are required to change the total number of particles in 
order to give rise to a chemical equilibrium. However these processes are suppressed 
by the VEV dependent masses of the gauge bosons (and gauginos). As a result the 
Universe enters a period of a quasi-thermal phase during which the reheat plasma 
evolves adiabatically. This lasts until the above mentioned scatterings become ef¬ 
ficient, at which time full equilibrium is finally attained, i.e., within one Hubble 
time. 

Such a late thermalization also results in a low reheat temperature. For the 
intermediate values of the flat direction’s VEV, the reheat temperature within MSSM 
varies: 

10 3 GeV < T R < 10 7 GeV. (1.4) 

Interestingly such a low reheat temperature does not lead to a significant generation 
of the gravitinos 9 . On the other hand, this is a bad news for thermal leptogenesis 


8 The flat directions could also reheat the Universe and could also explain the fluctuations present 
in the cosmic microwave background radiation, see |2l], |28). However there is a distinction, in this 
paper we demand that the inflaton decay is the major source for the entropy production. 

9 Note that gravitinos are still produced during preheating 0 |32[ and the final stage of 
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which typically requires that Tr > 10 J GeV [38 


One of the aims of this paper is to explore an alternative prescription for lepto- 
genesis, here, we will introduce a new paradigm, quasi-thermal leptogenesis, which 
points towards a successful baryogenesis within SUSY. 

The rest of this paper is organized as follows. In Section 2 we briefly review ther- 
malization after the perturbative inflaton decay (including relevant processes which 
lead to kinetic and chemical equilibrium), and then discuss the impact of supersym¬ 
metric flat directions on thermalization. Through numerical examples, in Section 3, 
we will underline the dramatically altered picture which emerges. In Section 4 we 
will consider thermalization in some of the other supersymmetric extensions of the 
SM, i.e., models with gauge mediation and split supersymmetry. After discussing 
particle production in the quasi-thermal phase in Section 5, we will specialize to 
the gravitino production and the leptogenesis in Sections 6 and 7 respectively. In 
Section 8 we will make some remarks on thermalization after preheating. We have 
also included Appendix which would help the paper to be self contained. 


2. Thermalization in supersymmetric theories: Part-1 

2.1 A brief recourse to thermalization 

For a plasma which is in full thermal equilibrium, the energy density, p, and the 
number density, n, of a relativistic particles are given by 

p = (tt 2 /30) T 4 , n — (C(3)/tt 2 ) T 3 , (Boson), 

p = (7/8) (tt 2 /30) T\ n= (3/4) (C(3)/tt 2 ) T 3 , (Fermion), (2.1) 

where T is the temperature of a thermal bath. Note that in a full equilibrium 
the relationships, ( E) ~ p 1 / 4 , and n ~ p 3 / 4 hold, with (E) = (p/n) ~ 3 T being 
the average particle energy. On the other hand, right after the inflaton decay has 
completed, the energy density of the Universe is given by: p ~ 3(FdMp) 2 . For a 
perturbative decay, which generates entropy, (see Eq. (|11.6|) of Appendix |11.1|) , we 
have ( E ) « m^ p 1 / 4 . Then, from the conservation of energy, the total number 

density is found to be, n & ( p/mfl) -C p 3 / 4 . Hence the complete inflaton decay results 
in a dilute plasma which contains a small number of very energetic particles. This 
implies that the Universe is far from full thermal equilibrium initially, 
perturbative decay of the inflaton jT| (see also |3(| [57j ). We will consider these cases separately 
in a different publication. 


7 










Reaching full equilibrium requires re-distribution of the energy among differ¬ 
ent particles, kinetic equilibrium , as well as increasing the total number of parti¬ 
cles, chemical equilibrium. Therefore both the number-conserving and the number- 
violating reactions must be involved. Full equilibrium is achieved shortly after the 
number-violating processes become efficient. From then on, we have the familiar 
hot big bang Universe, which consists of elementary particles in full thermal equilib¬ 
rium. The maximum temperature of the Universe after complete thermalization is 
the definition of the reheat temperature, i.e., Tr. 

Let us first begin with the number-conserving reactions which build kinetic equi¬ 
librium among the SM fermions. The most important processes are 2 —> 2 scatterings 
with gauge boson exchange in the f-channel, shown in Fig. (1). The cross-section for 
these scatterings is ~ a:|f| _1 . Here "t" is related to the exchanged energy, A E, and 
the momentum, Ap, through t = AU 2 —|Ap| . The fine structure constant is denoted 
by a (note that a > 10~ 2 in the MSSM). This cross section can be understood as 
follows: the gauge boson propagator introduces a factor of \t\~ 2 , while phase space 
integration results in an extra factor of \t\. Due to an infrared singularity, these 
scatterings are very efficient even in a dilute plasma 10 . 

Note that every degree of freedom of the MSSM has some gauge interactions 
with all other fields. Therefore any fermion in the plasma has t —channel scatterings 
off those fermions with the largest number density. For this reason the total number 
density, n, enters while estimating the scattering rate. 


10 There are also 2 —> 2 scattering diagrams with a fermion or scalar exchange in the t-channel. 
Diagrams with a fermion (for example, gaugino) exchange have an amplitude oc |t| , which will 

be canceled by the phase space factor |t|. This results in a much smaller cross-section oc s -1 , where 
s ~ 4 E 2 is the squared of the center-of-mass energy. Diagrams with a scalar (for example, Higgs) 
exchange are also suppressed by the following reason. A fermion-fermion-scalar vertex, which arises 
from a Yukawa coupling, flips the chirality of the scattered fermion. For relativistic fermions, as we 
consider here, the mass is <C E and a flip of chirality also implies a flip of helicity. This is forbidden 
by the conservation of angular momentum for forward scatterings, i.e., where t —» 0. As a result, 
the diagram in Fig. (2) has no t-channel singularity at all. Note that it is additionally suppressed 
by powers of Yukawa couplings compared to the diagram with gauge interactions in Fig. (1). For 
more details, see Ref. (|fj 



h 


h 


h 


~h 


Fig. 1: Typical scattering diagram which builds kinetic equilibrium in the reheat 
plasma. Note that the t —channel singularity which results in a cross-section oc |t| _1 . 


In addition one also needs to achieve chemical equilibrium by changing the num¬ 
ber of particles in the reheat plasma. The relative chemical equilibrium among dif¬ 
ferent degrees of freedom is built through 2^2 annihilation processes, occurring 
through s —channel diagrams. Hence they have a much smaller cross-section ~ as -1 . 
More importantly the total number of particles in the plasma must also change. It 
turns out from Eq. ( |2.1| ) that in order to reach full equilibrium, the total number of 
particles must increase by a factor of: n eq /n, where n & p/m$ and the equilibrium 
value is: n eq ~ p 3 / 4 . This can be a very large number, for examples given in the 
appendix, n eq /n ~ 0(1O 3 ). This requires that the number-violating reactions such 
as decays and inelastic scatterings must be efficient. 


Decays (which have been considered in Ref. pil) are helpful, but in general they 


cannot increase the number of particles to the required level. It was recognized in [^5| 
see also |29], 0, that the most relevant processes are 2^3 scatterings with gauge- 


boson exchange in the t— channel. Again the key issue is the infrared singularity of 
such diagrams shown in Fig. (2). The cross-section for emitting a gauge boson, whose 
energy is |f| i/2 <C E, from the scattering of two fermions is ~ a 3 |t| 1 . When these 
inelastic scatterings become efficient, i.e., their rate exceeds the Hubble expansion 


rate, the number of particles increases very rapidly [|tj, because the produced gauge 
bosons subsequently participate in similar 2 —> 3 scatterings. 


As a result the number of particles will reach its equilibrium value, n eq , soon 
after the 2^3 scatterings become efficient. At that point the scatterings and 
inverse scatterings occur at the same rate and the number of particles will not increase 
further. Therefore full thermal equilibrium will be established shortly after the 2 —> 3 
scatterings become efficient. For this reason, to a very good approximation, one can 
use the rate for inelastic scatterings as a thermalization rate of the Universe T t h r - 
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Fig. 2: Typical scattering diagrams which increase the number of particles. 

In order to estimate thermalization rate, we need to choose an infrared cut-off on 
the parameter t. With a reasonable choice of a cut-off, it turns out that in realistic 
models the 2 —> 3 scatterings have a rate higher than or same as that of the inflaton 


decay rate (for more details we refer the readers to |£5], ^9|, |30[). Therefore, if the 


inflaton decay products have gauge interactions, the Universe reaches full thermal 
equilibrium immediately after the inflaton decay. The reason is that the 2 —> 3 
scatterings with gauge boson exchange in the t —channel are very efficient. However 
this analysis crucially depends on having massless gauge bosons. 

In general the reheat plasma induces a mass for all particles which it contains. 
This is however negligible in a dilute plasma as in the case after the perturbative 
inflaton decay. However SUSY alters the situation quite dramatically as we shall 
discuss below. 


2.2 Flat directions in Supersymmetry 

The field space of supersymmetric theories contains many directions along which the 
D— and F —term contributions to the scalar potential identically vanish in the limit 
of unbroken SUSY. The most interesting such directions are those made up of SUSY 
partners of the SM fermions, namely the squarks and sleptons, and the Higgs fields. 
These directions have gauge and Yukawa interactions with matter fields, and hence 
decay before BBN. Therefore they do not lead to the cosmological moduli problem. 
As a matter of fact they can have very interesting cosmological consequences (for a 


review, see Ref. |26 


The superpotential for the Minimal Supersymmetric Standard Model (MSSM) 
is given by, see for instance pi 


kk mssm KQH u u T X d QH d d -1- \ e LH d e Y fiH u H d , 


( 2 , 2 ) 
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where H u , Hd, Q, L,u, d, e in Eq. (|2.2|) are chiral superfields representing the two 
Higgs fields (and their Higgsino partners), LH (s)quark doublets, RH up- and down- 
type (s)quarks, LH (s)lepton doublets and RH (s)leptons respectively. The dimen¬ 
sionless Yukawa couplings X u , A d, A e are 3x3 matrices in the flavor space, and we 
have omitted the gauge and flavor indices. The last term is the /i term, which is a 
supersymmetric version of the SM Higgs boson mass. 

The SUSY scalar potential V is the sum of the F- and D-terms and reads 

v = +\'£& D " D " (2.3) 

i a 


where 


Fi = 


dWi 


MSSM 


dxi 


D a = X *T° jXj ■ 


(2.4) 


Here the scalar fields, denoted by y, transform under a gauge group G with the 
generators of the Lie algebra and gauge coupling are by T a and g a respectively. 

For a general supersymmetric model with N chiral superfields X tl it is possible 
to find out the directions where the potential in Eq. (|2.3|) vanishes identically by 
solving simultaneously 


D a = X ] T a X = 0 , 


F = — = n 
dXi 


(2.5) 


Field configurations obeying Eq. (|2.5|) are called respectively D-flat and F-flat. 

D-flat directions are parameterized by gauge invariant monomials of the chiral 
superfields. A powerful tool for finding the flat directions has been developed in 
[T]. [T]]. g7| j where the correspondence between gauge invariance and flat 


directions has been employed. There are nearly 300 flat directions within MSSM 
The flat directions are lifted by soft SUSY breaking mass term, mo, 


V~m 2 0 \<p | 2 , 


as well as higher order terms arising from the superpotential [fh| 11 


( 2 . 6 ) 


2.3 Flat directions during and after inflation 

Because it does not cost anything in energy during inflation, where the Hubble 
expansion rate is Hi mo, quantum fluctuations are free to accumulate (in a 

11 Within supergravity there are corrections to the SUSY potential arising from the Kahler po¬ 
tential and also the mixing between superpotential and Kahler terms. We will consider these issues 
in Section [|. 


11 








coherent state) along a flat direction and form a condensate with a large VEV, po- 
Because inflation smoothes out all gradients, only the homogeneous condensate mode 
survives. However, the zero point fluctuations of the condensate impart a small, and 
in inflationary models a calculable, spectrum of perturbations on the condensate [^] . 
In an abuse of language we will collectively call such condensates as flat directions. 

After inflation, H oc £ -1 , the flat direction stays at a relatively larger VEV due to 
large Hubble friction term, note that the Hubble expansion rate gradually decreases 
but it is still large compared to mo. When H ~ mo, the condensate along the flat 
direction starts oscillating around the origin with an initial amplitude ~ po- From 
then on \p\ is redshifted by the Hubble expansion oc H for matter dominated and 
~ i/ 3 / 4 f or radiation dominated Universe. 

If higher order superpotential terms are forbidden, due to an R —symmetry (or a 
set of R— symmetries) |48|], then we naturally have, po ~ Mp |45|]. On the other hand, 
p 0 -C Mp will be possible if non-renormalizable superpotential terms are allowed. 
Let us first analyze the model independent case where we treat p as a free parameter 
which can vary in a wide range with an upper limit p o < Mp, and study various 
consequences. Note that a lower bound Hi < p 0 is set by the uncertainty due to 
quantum fluctuations of p during inflation. 


2.4 Flat directions and inflaton decay 


If a flat direction which has a non-zero VEV has couplings to the inflaton decay 
product(s), then it will induce a mass, y\p\, where y is a gauge or Yukawa coupling. 
The inflaton decay at the leading order will be kinematically forbidden if y\p\ > 
m^/2. One should then wait until the Hubble expansion has redshifted, \p\, down 
to (m r j ) /2y). The decay happens when (note that \p\ oc H , after the flat direction 
starts oscillating and before the inflaton decays): 


Hi = min 




m 0 ,T d 


(2.7) 


LVWo, 

The inflaton also decays at higher orders of perturbation theory to particles which 
are not directly coupled to it 0- This mode is kinematically allowed at all times, 

m^/y |v?|) 2 T d . It will become efficient at 

2/3 


but the rate is suppressed by a factor of 

f mpmo\ i/a 

2 ~ l ) d ' 

Therefore, if the decay products are coupled to a flat direction with a non-zero VEV, 
the inflaton will actually decay at a time when the expansion rate of the Universe is 


( 2 . 8 ) 
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given by 


H d = max [Hi, H 2 ] . (2.9) 

In general it is possible to have, H d <C Td, particularly for large values of <p o 12 . 
Flat directions can therefore significantly delay inflaton decay on purely kinematical 
grounds 13 . 

2.5 Flat directions and thermalization 

Flat directions can dramatically affect thermal history of the Universe even if they do 
not delay the inflaton decay. The reason is that the flat direction VEV spontaneously 
breaks the SM gauge group. The gauge fields of the broken symmetries then acquire a 
supersymmetry conserving mass, m g ~ g\tp\, from their coupling to the flat direction, 
where g is a gauge coupling constant. 

The simplest example is the flat direction corresponding to the H u H d monomial. 
One can always rotate the field configuration to a basis where H n X = H d2 = 0, 
with subscripts 1 and 2 denoting the upper and lower components of the Higgs 
doublets respectively. The complex scalar field, p = (. H u 2 + H d i) /a/2, represents a 
flat direction. It breaks the SU(2)w x U(l)y down to f/(l) em (exactly in a similar 
fashion as what happens in the electroweak vacuum). The W ± and Z gauge bosons 
then obtain a mass from their couplings to the Higgs fields via covariant derivatives. 

The complex scalar field, (H u j + H d2 ) /\/2, and the real part of, (H Uj2 — H dt i) / a/2, 
also acquire the same mass as W ± and Z, respectively, through the D— term part of 
the scalar potential. The Higgsino fields H u X and H d 2 are paired with the Winos, 
while (ti U)2 — H d ^j /y/2 is paired with the Zino, and acquire the same mass as W ± 
and Z, respectively, through the gaugino-gauge-Higgsino interaction terms. In the 
supersymmetric limit, the flat direction and its fermionic partner (^H Ut2 + H dj ij j\[2, 
as well as the photon and photino, remain massless. 

Many of the flat directions break the entire SM gauge group. The prominent ex¬ 
amples are flat directions corresponding to the LLddd and QuQue monomials. The 
whole SM gauge group can be also broken by multiple flat directions which individ- 
12 Flat directions can also affect inflaton decay in other ways (4tJ. 

13 Note that such a delayed decay can also naturally implement the modulated fluctuations mech¬ 
anism for generating adiabatic density perturbations by converting their isocurvature perturba¬ 
tions |b|]. Due to the dependence of H d on <p 0 , fluctuations of the flat direction (~ Hi) imply 
an inhomogeneous inflaton decay which can give rise to perturbations of the correct magnitude if, 
Hi ~ 10"Vo- 
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ually break only part of the symmetry. For example flat directions corresponding to 
the LLe monomial 14 break the SU(2)w x U(l)y group, while directions correspond¬ 
ing to the udd monomial break the SU(3)c- Note that all independent flat directions 
can simultaneously acquire a large VEV during inflation 15 . 

Let us now imagine one of such flat directions which has obtained a large VEV 
during inflation. In which case a flat direction can crucially alter thermal history 
of the Universe by suppressing thermalization rate of the reheat plasma. Note that, 
m g , provides a physical infrared cut-off for scattering diagrams with gauge boson 
exchange in the t —cha nn el shown in Figs. (1, 2). The thermalization rate will then 
be given by (up to a logarithmic “bremsstrahlung” factor): 

Tl 

Tthr ~ a 2 -— ~ 2 , (2.10) 

m 

where we have used m 2 ~ ct|<y3| 2 . After the flat direction starts its oscillations at 
H ~ mo, the Hubble expansion redshifts, \tp\ 2 oc R~ 3 , where R is the scale factor of 
the FRW Universe. The interesting point is that, n oc R~ 3 , as well, and hence T t hr 
remains constant while H decreases for H < mo- This implies that r t h r eventually 
catches up with the expansion rate, even if it is initially much smaller, and shortly 
after that the full thermal equilibrium will be achieved. 

Therefore the flat directions, even if they do not delay the inflaton decay, will 
modify thermal history of the Universe. Depending on whether mo > Td or mo < 
Td, different situations will arise which we discuss separately. If two or more flat 
directions with non-zero VEVs induce mass to the gauge bosons, then \ip\ denotes 
the largest VEV. 


• m 0 > r d : 

In this case the inflaton decays after the flat direction oscillations start. The 

14 If the LLe monomial develops a VEV during inflation, which breaks the electroweak symmetry, 
it generates a mass to the weak hypercharge gauge boson which mixes to the photon. It is then 
possible to excite that gauge boson from vacuum fluctuations, which would then be converted to 
the primordial magnetic field, see |52| ]. 

15 Only exceptional flat direction is the H u Hd direction. Since it has superpotential Yukawa 
couplings to all MSSM fields, it will induce a large mass to other scalar fields. As a result no other 
flat direction can develop a large VEV in the presence of H u Hd . This is not the case for the H U L 
flat direction: LLddd has no Yukawa couplings to it, and hence can simultaneously obtain a large 
VEV. For a discussion on multiple flat directions, see [j53f. 
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inflaton oscillations, which give rise to the equations of state close to non- 
relativistic matter, dominate the energy density of the Universe for H > T d . 
This implies that R oc H~ 2 / 3 , and \(p\ is redshifted oc H in this period. We 
therefore have, 


<,Pd ~ 



( 2 . 11 ) 


where (p$ denotes the amplitude of the flat direction oscillations at the time 
of the inflaton decay H ~ T d . The number density of particles in the reheat 
plasma at this time is given by: 


iot|m| 

rid - 

m# 


( 2 . 12 ) 


where we have used p ~ 3 (T d Mp) 2 . Note that \ip\ and n are both redshifted 
oc R 3 for H < T d . Thus, after using Eq. ( p,10|) , we find that complete 
thermalization occurs when the Hubble expansion rate is 

2 


H< 


thr 


10 a 2 


Mv 


m, 


o 


V o J rn ^ 


(2.13) 


• m 0 < T d : 

In this case the flat direction starts oscillating after the completion of inflaton 
decay. The Universe is dominated by the relativistic inflaton decay products 
for H < T d , implying that R oc H. 

The number density of particles in the plasma is redshifted oc H 3 / 2 and, see 
Eq. (|T1), at H = mo, it is given by 


10 T\M 2 

n 0 - 

m ^ 



(2.14) 


Note that n, \ip\ 2 (x R 3 , and hence T t hr remains constant, for H < mo- The 
reheat plasma then fully thermalizes when the Hubble expansion rate is 

1 f M P \ 2 ml 




thr 


lOct 


m 0 


V J 


(2.15) 


Note that the kinetic equilibrium is built through 2 —> 2 scattering diagrams 
as in Fig. (1), which have one interaction vertex less than those in Fig. (2). This 
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implies that the rate for establishment of kinetic equilibrium will be I\i n ~ ct 1 r t h r - 
Therefore we have a typical relationship: 


r c i r kin > r 


thr 


(2.16) 


in SUSY 16 . 

This implies that the Universe enters a long period of quasi-adiabatic evolution 
after the inflaton decay has completed. During this phase, the comoving number 
density and (average) energy of particles remain constant. It is only much later 
that the 2 —> 2 and 2 —> 3 scatterings become efficient and the plasma completely 
thermalizes. As we shall see, this has important cosmological consequences. 

One comment is in order before closing this section. So far we have neglected 
the decay of the flat directions and their interactions with the reheat plasma. These 
effects are considered in an Appendix |11.4 , and shown to be negligible before the 
Universe fully thermalizes. 


3. Reheat temperature of the Universe 


The temperature of the Universe after it reaches full thermal equilibrium is referred 
to as the reheat temperature Tr. In the case of SUSY, we therefore have: 

T r ~ (H thl M P ) 1/2 , (3.1) 


where, depending on the details, H thr is given by Eqs. (|2.13|) and Q2.15| ). 

As mentioned earlier, we typically have H thv <C T d . Therefore in SUSY theo¬ 
ries the reheat temperature is generically much smaller than the standard expression 
Tr ~ (TdMp) 1 / 2 , which is often used in the literature and assumes immediate ther- 
malization after the inflaton decay. 

An interesting point to note that the reheat temperature depends very weakly 
on the inflaton decay rate, for instance Eq. (|2.15| ) implies that Tr oc T^/ 4 , while 
Tr is totally independent of T d in Eq. (|2.13|) . This is not difficult to understand, 
regardless of how fast the inflaton decays, the Universe will not thermalize until the 
2 —> 3 scatterings become efficient. The rate for these scatterings essentially depends 
on the flat direction VEV and mass. 


16 Relative chemical equilibrium among different degrees of freedom is built through 2 —> 2 anni¬ 
hilations in the s—channel with a rate ~ a 2 n/E 2 -C r t hr- Hence composition of the reheat plasma 
will not change until full thermal equilibrium is achieved. 
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As expected, a larger p$ results in slower therma li zation and a lower reheat 
temperature. On the other hand, larger values of m 0 lead to a higher Tr. Since 
the flat direction oscillations start earlier in this case, its VEV (thus mass of gauge 
bosons) is redshifted faster and thermalization rate will be less suppressed. 

The reheat temperature, Tr, depends on the inflation sector through the inflaton 
mass m<y It is counterintuitive, but the fact is that, from the above Eqs. Q2.131 ), ( P ■ 1 5| ) 
and 0 , a larger m^ results in a lower reheat temperature. This is due to the con¬ 
servation of energy which implies that the number density of inflaton decay products 
is inversely proportional to their energy, which is ~ initially. Since the scattering 
rate is proportional to the number density, a larger thus results in a smaller T t h r , 
and a lower Tr. 

If Hthv is very large or T d is very small, we may find H thr > T d . Obviously 
thermalization cannot occur before the inflaton decay has completed. This merely 
reflects the fact that the 2 —> 3 scatterings are already efficient when H ~ T d . This 
will be the case if the flat direction VEV is sufficiently small at the time of the 
inflaton decay, and/or if the reheat plasma is not very dilute. The former happens 
for a small <po or large mo, while smaller values of m^ lead to the latter case. The 
reheat temperature in such cases follows the standard expression: Tr ~ (T&M-p) 1 ^ 2 . 


3.1 Numerical examples 


We now present some examples to demonstrate the impact of flat directions on the 
rate of thermalization. We choose the nominal value of the inflaton mass, m^ = 
10 13 GeV, typical value for the flat direction mass in models with a weak scale SUSY, 
m 0 — 100 GeV —1 TeV, and four representative VEVs, p 0 = M P , 10~ 2 M P , 10 _4 M P , 
and po < 10 -6 M P . 

If the inflaton decays gravitationally, see the example in Appendix |11.2| , we have 
T d ~ 10 GeV. In this case H t h r is given by Eq. (|2.13| ). For larger inflaton couplings 
to matter, T d > mo can be obtained, in which case H t h r follows Eq. (|2.15|) . As a 
sample case we have chosen T d = 10 4 GeV. 
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VEV (in GeV) 

T R (T d = 10 GeV) 

T R (T d = 10 4 GeV) 

<Po ~ Mp 

3 x 10 3 

7 x 10 4 

ipo = 10 2 Mp 

3 x 10 5 

7 x 10 6 

ipo = 10 4 Mp 

3 x 10 7 

7 x 10 8 

( p 0 < 10 -6 Mp 

3 x 10 9 

7 x 10 10 


Table 1: The reheat temperature of the Universe for the inflaton mass, = 10 13 
GeV, and two values of the inflaton decay rate, T d = 10, 10 4 GeV. The flat direction 
mass is m 0 ~ 1 TeV. The rows show the values of Tj> for flat direction VEVs varying 
in a wide range. 

It is evident that thermalization rate (hence Tp) becomes smaller as tpo increases. 
It is remarkable that in the extreme case where <^ 0 ~ Mp, the reheat temperature can 
be as low as TeV. This is in stark contrast with the case if thermalization was instant 
(as expected in a non-SUSY case), which would result in hierarchically higher reheat 
temperatures Tp > 10 9 GeV (for the chosen values of r d ). Note that thermalization 
still remains quite slow for ipo rsj 10 4 Mp. It is only for ip o < 10 6 Mp that the flat 
direction VEV is sufficiently small in order not to affect thermalization, thus leading 
to the standard expression Tp ~ {Y^Mp) 1 ^ 2 . 

This clearly underlines the fact that complete thermalization can be substantially 
delayed in supersymmetry. Indeed, within the range determined by the uncertainty 
of the quantum fluctuations, < p 0 < M P , flat directions considerably slow down 
thermalization. This has important cosmological consequences which will be consid¬ 
ered in detail in the following sections. 


4. Thermalization in supersymmetric theories: Part-2 

So far we have not made any specific assumption on the mediation of SUSY breaking 
to the observable sector, or the nature of higher order terms which lift the flat 
directions. Instead we considered the typical case in models with soft masses at TeV 
scale, and we treated the flat direction VEV as a free parameter which is bounded 
from above by Mp. In this section we shall study issues and some subtleties that 
may arise in more detail. 




4.1 Higher order superpotential terms and Kahler corrections 


In models with gravity [|4]] and anomaly |>4] mediation SUSY breaking results in as 
usual soft term, mg| 9 ?| 2 , in the scalar potential where mo — 100 GeV — 1 TeV. There 
is also a new contribution arising from integrating out heavy modes beyond the scale 
M, which usually induces non-renormalizable superpotential term: 


W ~ A,- 


<j> r 


M n -s ’ 


(4.1) 


where 4> denotes the superfield which comprises the flat direction p. In general M 
could be the string scale, below which we can trust the effective field theory, or 
M = Mp (in the case of supergravity). In addition, there are also inflaton-induced 
supergravity corrections to the flat direction potential. By inspecting the scalar 
potential in N = 1 supergravity, one finds the following terms 




(k^k^k^ 


A4p J ’ 


(k^K'^Dj 


W*(I) 

Ml 



1 


(4.2) 


where Dj = d/dl + KjW/M'l. The Kahler potential for the flat direction and the 
inflaton are given by K(p*,p) and K(I*,I), and W(I) denotes the superpotential 
for the inflaton sector. 

All these terms provide a general contribution to the flat direction potential 
which is of the form ||45|| : 


= H'tHf ( ^ 


<p 


(4.3) 


where / is some function. Such a contribution usually gives rise to a Hubble induced 
correction to the mass of the flat direction with an unknown coefficient, which de¬ 
pends on the nature of the Kahler potential. The relevant part of the scalar potential 
is then given by |k) 


UD ( m 2 0 + c H H 2 ) \p\ 2 + A, 


\<P\ 


2(n—1) 


M 2 ( n -3) ’ 


(4.4) 


with n > 4. Note that c# can have either sign. If > 1, the flat direction mass 
is > H. It therefore settles at the origin during inflation and remains there. Since 
| ip | = 0 at all times, the flat direction will have no interesting consequences in this 
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case 17 . The case with ch < 0 will be more interesting. A negative ch can arise at a 
tree-level [f|5[ , or as a result of radiative correction |]56|, [57| . The flat direction VEV 


is in this case driven away from the origin and quickly settles at a large value which 
is determined by higher order terms that stabilize the potential. A large VEV can 
also be obtained if 0 < ch < 1 . 

It is possible to eliminate certain or all higher order superpotential terms terms 
by assigning a suitable R —symmetry (or a set of A—symmetries) |48|]. However let us 
assume that all such terms which respect the SM gauge symmetry are indeed present. 


As shown in [[47[, all of the MSSM flat directions are lifted by higher-order terms 
with n <9. If a flat direction is lifted at the superpotential level n, the VEV that 
it acquires during inflation will be: tpi rv./ (HjM n 3 ) 1 / n 2 i where Hi is the expansion 
rate of the Universe in the inflationary epoch. If ch < 0, this will be the location of 
the minimum of the potential stabilized by the higher-order term. 

After inflation, the flat direction VEV slides down to an instantaneous value: 
(H(t)M n ~ 3 ) 1 ^ n 2 . Once H(t ) ~ m 0l soft SUSY breaking mass term in the potential 
takes over and the flat direction starts oscillating around its origin with an initial 
amplitude: <p 0 ~ (moM" -3 ) 1 ^ 2 . If M = Mp, we will have ipo ~ 10 10 GeV for n = 4, 
and (fi 0 ~ 10 16 GeV for n — 9. In particular, tp 0 > 10 14 GeV, for the flat directions 
lifted at n > 6 levels. These directions affect thermalization considerably and lower 
the reheat temperature of the Universe as we discussed earlier, see Table. 1. 


4.2 Models with gauge-mediated supersymmetry breaking 


Our estimations of the thermalization time scale, Eqs. Q2.13 ) and (|2.15|) , are valid if 
the flat direction potential is quadratic ml\ip\ 2 for all held values (up to Mp). However 


the situation is more subtle in models with gauge-mediated SUSY breaking |^8] 


In these models there exists a sector with gauge interactions that become strong 
at a scale A DS b -C Mp. This induces a non-perturbative superpotential and leads to a 


17 The positive Hubble induced mass to the flat direction has a common origin to the Hubble in¬ 
duced mass correction to the inflaton in supergravity models. This is the well known 77 -problem (55), 
which arises because of the canonical form of the inflaton part of the Kahler potential. A large 
77 generically spoils slow roll inflation. In order to have a successful slow roll inflation, one needs 
77 <C 1 (and hence ch -C 1). Note that the origin of the Hubble induced corrections to the mass of 
flat direction is again the cross terms between inflaton Kahler potential and superpotential terms 
and flat direction Kahler potential. There is no absolutely satisfactory solution to the 77 problem 
and therefore large coefficient, i.e., ch ~ +0(1). In this paper we assume that somehow the 77 
problem has been addressed and therefore also ch +0(1). 
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non-zero F-coniponent, ~ A BSB , for a chiral superfield which breaks supersymmetry. 
The gravitino mass will therefore be: 7713/2 rv./ A BSB /M P . At a next step, SUSY 
breaking is fed to a messenger sector with a mass scale m mess . Finally, the soft SUSY 
breaking parameters are induced in the observable sector by integrating out the 
messenger fields, which have some common gauge interactions with the observable 
sector, resulting in soft scalar masses m 0 . 

As a result, the flat direction potential has the following forms: 

V = m 2 0 \<p\ 2 

V ~ (m 0 ?77mess) 2 


In 


M 


m v 


V = ml/M' 


^mess Yl \v\ Y 

M > 


|/91 Y 777 mess , 

mQTn mess 


m 3/ 2 
777o777 m ess 
m 3 / 2 


(4.5) 


This behavior can be understood as follows. The flat direction obtains its soft SUSY 
breaking mass by integrating out the messenger fields. For \tp\ > M mess these fields 
acquire a large mass oc \tp\ from the flat direction VEV, for details, see 0- The 
flat direction mass, i.e., y/V", then decreases and the potential varies very slowly. 
At very large field values, however, dominance of the contribution from the gravity 
mediation recovers the conventional, |/9| 2 , dependence on the potential, although now 
the gravity mediated soft SUSY breaking contribution is 777,3/2 instead m 0 . 

Now let us find the thermalization time scale in models with gauge mediation. 
For small values of 7773 / 2 , typically arising in these models, the flat direction starts 
its oscillations when H = m 3 / 2 and U(/9 0 ) — m l/ 2 L Po- Also, the oscillations usually 
start after the inflaton has decayed, i.e., such that T d > 7773 / 2 . 

Note that \(f\ 2 and the number density, 77, are redshifted alike oc R~ 3 for H < 
7773/2, and hence T t hr remains constant, while H decreases. This changes when |</?| < 
m o m mess/ ? » 3 / 2 ) after that \ip\ is redshifted oc R~ 3 ||2TJ, and hence r thr increases oc i? 3 . 
The transition occurs when H ~ 777 0 777 mess //9 0 - Here we have assumed that the 
Universe is radiation dominated, thus R oc H ~ 1 ^ 2 , for H < r d . Due to a rapid 
increase of the ratio: r t h r/H, the reheat plasma reaches full thermal equilibrium 
when |/91 > 777 mess . After using Eqs. ( p.lC| ) and (|2.14|) , we find that the Hubble 
expansion rate at the time of thermalization is then given by: 
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thr 
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(4.6) 


The maximum impact on thermalization happens in models with low-energy gauge 
mediation [ETJ, which can give rise to very light gravitinos. In these models A B g B , 


21 
















rsj 


m mess and m 0 are all related to each other by one-loop factors: m mess ~ Adsb/16vt 2 
and mo ~ ^mess/16vr 2 [^l|. After using these relations and taking into account that 
m 3/2 ~ Aq SB /M p , Eq. ( TO reads 

1/5 




thr 
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r d m 3/ 2 \ (Ml 

m l J 


<Po 


7/5 


™3/2- 


(4.7) 


The reheat temperature then follows: Tr ~ (H^Mp) 1 ^ 2 . 

For example, consider a model where A DSB ~ 10' GeV, m mess ~ 10 5 GeV and 
mo ~ 1 TeV. The gravitino mass in this model is m. 3/2 — 100 KeV. Note that 
the potential is quadratic for \ip\ < 10 3 GeV and \ip\ > 10 12 GeV, while being 
logarithmic in the intermediate region. Table. 2 summarizes the reheat temperature 
in this model for similar values of m^, T d and </?o as in Table. 1. What we find is 
that T b is considerably lower in this case. 


VEV (in GeV) 

T R (T d = 10 GeV) 

T R (T d = 10 4 GeV) 

(po rv - / Mp 

8 x 10 2 

3 x 10 3 

ipo = 10 - 2 M P 

2 x 10 4 

6 x 10 4 

ipo = 10 ” 4 M P 

5 x 10 5 

2 x 10 5 

ipo = 10 ^ Mp 

10 7 

3 x 10 7 


Table 2: The reheat temperature of the Universe in a model with low-energy gauge 
mediation where the gravitino mass is m.3/2 = 100 KeV. The inflaton mass and its 
decay rate are the same as in Table. 1. The rows show the values of Tr for flat 
direction VEVs chosen the same as in Table. 1. 


To conclude, in models with gauge mediation the Universe thermalizes even 
more slowly and the resulting reheat temperatures are even lower compared to that 
of gravity mediation case. This is mainly due to the fact that the flat directions start 
their oscillations much later in this case. 


4.3 Split supersymmetry 

Now we consider the opposite situation where mo 3> 1 TeV. This happens in the re¬ 


cently proposed split SUSY scenario [Q. This scenario does not attempt to address 
the hierarchy problem. It allows the scalars (except the SM Higgs doublet) to be 
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very heavy, while keeping the gaugino and Higgsino fields light. This removes prob¬ 
lems with flavor changing and CP —violating effects induced by scalars at one-loop 
level. On the other hand, it preserves attractive features like supersymmetric gauge 
coupling unification and a light LSP which can account for the dark matter. 


Note that from Eqs. (|2.13|) and ( |2.15|) , the thermalization time scale becomes 
shorter for larger values of mo, because when mo is larger the flat directions start 
oscillating earlier. Their VEV and the induced masses for the gauge bosons will 
therefore be redshifted faster in this case. Let us consider the favored range of the 
soft scalar masses in split SUSY: 10 8 GeV < mo Y 10 13 GeV. This removes the flavor 
changing and CP —violating effects and results in an acceptable gluino lifetime [6Cj 


For the inflaton mass rri^ = 10 13 GeV, we find H t h r > 100 GeV for all values of 
<Po < M P . 


Therefore in the case of split SUSY, we always have 


T r > 10 10 GeV. 


(4.8) 


This is fairly a robust prediction. 


This might seem unacceptably high for thermal gravitino production. However, 
models of split SUSY can also accommodate gravitinos with a mass m.3/2 > 50 
TeV jOH. |0T||. Such superheavy gravitinos decay before BBN fllTj, and hence are 


not subject to the tight bounds coming from BBN |19], |2(J. As a matter of fact, 
the gravitino overproduction can in this case turn into a virtue. The late decay 
of gravitinos, below the LSP freeze-out temperature, can produce the correct dark 
matter abundance in a non-thermal fashion f 
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Moreover, if produced abundantly, the gravitinos will dominate the Universe. 
Their decay in this case dilutes the baryon asymmetry. Consequently, larger parts of 
the parameter space will become available for thermal leptogenesis and flat direction 
baryogenesis p2l. 


We conclude that, due to larger scalar masses, thermalization is not affected by 
the flat directions in split SUSY. However, for m.3/2 > 50 TeV, having high reheat 
temperatures and thermal overproduction of gravitinos can be a privilege rather than 
a handicap |3^, . 
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5. Particle production during thermalization 


5.1 Quasi-adiabatic evolution of the Universe 


Right after the inflaton decay has completed the energy density of the Universe 
is given by, p ~ 3(TdMp) 2 , and the average energy of particles is (E) ~ m ( f > 18 . 
Deviation from full equilibrium can be quantified by the parameter "A" Jl|], where 


^ ~ j'i 


10 4 


TaMp 


m. 


(5.1) 


Here we dehne T ps {E)/ 3, in accordance with full equilibrium. Note that in full 
equilibrium, see Eq. o, we have A ~ g* (= 228.75 in the MSSM). On the other 
hand, see Eq. (|11.6|) , after the inflaton decay we have Td -C m|/Mp which implies 
that A -C 228.75. One can also associate parameter Ai = Spi/T* to the i —th degree 
of freedom with the energy density pi (all particles have the same energy E, and 
hence T, as they are produced in one-particle decay of the inflaton). Note that 
A = YhiAi, and in full equilibrium we have Aj — 1. 

As we have discussed in Section |2j. thermalization is very slow in supersymmetry. 
As a result, the Universe enters a long phase of quasi-adiabatic evolution during which 
the comoving number density and comoving (average) energy of particles in the 
plasma remain constant. Since particles are produced in one-particle inflaton decay, 
the distribution is peaked around the average energy. The 2-^2 scatterings (which 
become efficient shortly before complete thermalization) smooth out the distribution. 
This implies that the reheat plasma is in a quasi-thermal state which is not far from 
kinetic equilibrium, but grossly deviated from chemical equilibrium. In this period, 
which lasts until the 2 —> 3 scatterings shown in Fig. (2) become efficient, the Hubble 
expansion redshifts pt oc A” 4 , n* oc R~ 3 and T oc iU 1 . Therefore A and A t remain 
constant throughout this epoch. Note that A depends on the total decay rate of 
the inflaton Tj and its mass through Eq. O. While, Ai are determined by 
the branching ratio for the inflaton decay to the i —th degree of freedom 19 . The 
composition of the reheat plasma is therefore model-dependent before its complete 
thermalization. However, some general statements can be made based on symmetry 
arguments: 

ls For example, in a two-body decay of the inflaton, we have exactly E = m<p/2. 

19 We would like to thank Antonio Masiero for highlighting this fact. 
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• For CP-conserving couplings, the inflaton decay produces the same number of 
particles and anti-particles associated to a given held. 

• For a singlet inflaton (which is the case in almost all models) 20 the gauge 
invariance implies that inflaton has equal couplings (thus branching ratios) to 
the fields which are in an irreducible representation of the gauge (sub)group. 
This holds in the presence of spontaneous symmetry breaking via hat direction 
VEVs, as long as the particle masses induced by coupling to the hat direction 
are (much) smaller than that of the inflaton mass. 


• If the inflaton mass is (much) larger than the soft SUSY breaking masses in 
the observable sector, the inflaton decay produces the same number of bosonic 
and fermionic components of matter superhelds 21 . 

During the quasi-adiabatic evolution of the reheated plasma, i.e., for < H < 
r d , we have 22 


pi = Ai—T 4 , m = Ai—T A , 


71 


71 


and the Hubble expansion rate follows 

H ~ A 1/2 


T 2 


3 A/p 


(5.2) 


(5.3) 


In this epoch T varies in a range T m \ n < T < T max , where T max ~ m^/3 is reached 
right after the inflaton decay. Because of complete thermalization, T sharply drops 
from T m i n to Tr at H t h r , where the conservation of energy implies that 


Tr = 


A 


228.75 


1/4 


T 

± rr 


(5.4) 


20 Only exceptional model is an example of assisted inflation 163 driven by N supersymmetric 
flat directions Q which are not gauge singlets. 

21 This holds for a perturbative decay which is relevant for the last stage of reheating. Non- 
perturbative inflaton decay via preheating, which may happen at an initial stage, produces bosons 
much more abundantly than fermions as a result of strong SUSY breaking by large occupation 
numbers. 

22 0ne can express A, in terms of a negative chemical potential pi, where At = exp (pi/T). 
Note that for a large negative chemical potential, i.e., in a dilute plasma, the Bose and Fermi 
distributions are reduced to the Maxwell-Boltzmann distribution and give essentially the same 
result. The assignment of a chemical potential merely reflects the fact that the number of particles 
remains constant until the number-violating reactions become efficient. It does not appear as a 
result of a conserved quantity (such as baryon number) which is due to some symmetry. Indeed, 
assuming that inflaton decay does not break such symmetries, here we assign the same chemical 
potential to particles and anti-particles. 
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The (final) entropy density is given by s = (27t 2 /45) g*T^, where we take g * = 228.75. 


5.2 Quasi-thermal particle production 


An important cosmological phenomenon is the production of (un)stable particles 
arising in theories beyond the SM in the early Universe. Let us consider such a 
particle, y, with a mass, m x , which is weakly coupled to the (MS)SM fields. Through 
its couplings, y will be inevitably produced in a plasma consisting of only the (MS)SM 
particles. In full equilibrium, one can calculate the abundance of a given particle by 
taking a thermal average of the relevant processes (for example, see Ref. |B5|). The 
only inputs required are the particle mass and its couplings to the (MS)SM fields, 
plus the reheat temperature Tfi. This leads to a lower bound on the relic abundance 
which is independent from the details of the reheating (except for Tr). 

Here we are interested in particle production during the quasi-thermal phase of 
the Universe. Usually, scatterings are the most important processes for production 
of x■ The number density of y, denoted by n x , then obeys 


n. 


3Hn x = 


X^’ rel 

i<j 


(T; 




)niUj 


(5.5) 


Here n * and nj are the number densities of the i— th and j— th particles, (Jij^ x * s 
the cross-section for producing y from scatterings of i and j, and the sum is taken 
over all fields which participate in y production. Also () denotes averaging over 
the distribution. Note that production will be Boltzmann suppressed if T < m x / 3. 
Therefore, to obtain the total number of y quanta produced from scatterings, it will 
be sufficient to integrate the RH side of (|5.5| ) from the highest temperature down to 
m x j 3. The physical number density of y produced at a time t\ will be redshifted 
by a factor of (R 2 /R 1 ) 3 at any later time t 2 . If the Universe is filled by relativistic 
particles, H = (1/2 1) and R oc H -1 / 2 . The relic abundance of y, normalized by the 
entropy density, s, will then be 


!!* „ 10 —5 
s \ A 




A-iAj (t^rel ®ij —>x) -^P 


(5.6) 


where we have used Eq. (|R4l). 


6. Gravitino production in a quasi-thermal case 

As we discussed in Section 3, slow thermalization in SUSY generically results in 
a low reheat temperature compatible with the BBN bounds on thermal gravitino 
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production. However gravitinos are also produced during the quasi-thermal phase 
prior to a complete thermalization of the reheat plasma. There are various chan¬ 
nels for the gravitino production from scatterings of gauge, gaugino, fermion and 
sfermion quanta. The scattering cross-section for all such processes is oc 1/Mp, but 
the numerical coefficient depends on a specific channel. In full thermal equilibrium, 
all degrees of freedom (except for their bosonic or fermionic nature) have the same 
occupation number, and hence obtaining the total production cross-section is rather 
easy. However, as mentioned earlier, the composition of the reheat plasma is model- 
dependent in a quasi-thermal phase. Therefore calculating the total cross-section is 
more involved in this case. 

An important distinction is that we only need to consider scatterings of fermions 
and/or sfcrmions for the following reasons; the gauge and gaugino quanta have large 
masses ~ a 1//2 (^d (induced by the flat direction VEV) at a time most relevant for the 
gravitino production, i.e., when H ~ T^, therefore, they decay to lighter fermions and 
sfcrmions at a rate ~ Here a 3 ^ 2 p>d is the decay width at the rest frame 

of gauge/gaugino quanta, and Lpd/m^ is the time-dilation factor. The decay rate is 
Td, thus gauge and gaugino quanta decay almost instantly upon production, and 
they will not participate in the gravitino production 23 . This is indeed welcome as 
scatterings which include gauge and/or gaugino quanta in the initial state (particu¬ 
larly scattering of two gluons) have the largest production cross-section Jl4|, [U|. 

As a consequence, production of the helicity ±1/2 states will not be enhanced 
in a quasi-thermal phase as scatterings with a gauge-gaugino-gravitino vertex will be 
absent. 

Now let us discuss the relevant scattering processes and production cross-sections 
in the quasi-thermal phase 24 . The total cross-section, and cross-sections for mul- 
tiplets comprising the LH (s)quarks Q, RH up-type (s)quarks u , RH down-type 
(s)quarks d, LH (s)lcptons L , RH (s)leptons e and the two Higgs/Higgsino doublets 
H u , Hd have been given separately for each channel 25 . Here 1 < i.jj < 3, a, b = 1,2 
and 1 < a,/3 < 3 are the flavor, weak-isospin and color indices of scattering degrees 

23 The inverse decay happens at a rate ~ an/E 2 and, as mentioned earlier, it is inefficient before 
complete thermalization of the plasma. This is just a reflection of the fact that the reheat plasma 
is out of chemical equilibrium for H > Hthr- 

24 0ne might wonder that scatterings of (s)fermions off the flat direction condensate would be the 
dominant source as the latter carries a much larger number of (zero-mode) quanta. However, these 
scatterings are suppressed for the same reasons as discussed in Appendix |ll.4 . 

25 We also note that these cross-sections are free from infrared logarithmic divergences |lj, |l5||. 
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of freedom respectively. Also a 3 , a 2 and ay are the gauge fine structure constants 
associated to the SU(3)c, SU(2)w and U(l)y groups respectively. For each channel 
we have taken the total cross-section from the table on page [181] of the first refer- 
ra, and split it according to the gauge interactions for the relevant flavors 


ence m 


(and each of the two doublets in the case of Higgs). Note that the prefactor 3/4, 
which arises as a result of Fermi-Dirac statistics in full thermal equilibrium will not 
be relevant here. 


• fermion + anti-sfermion —» gravitino + gauge field ,, 

sfermion + anti-fermion —>■ gravitino + gauge field: The total cross-section for 
this channel is (1/32M|) x (48ay + 21a 2 + lloy), see row (E) of the table on 
page [181] of the first reference in [[15]. For different degrees of freedom we then 
find: 


Vi — a 2 — 


32 Mi SijX 


2 . 3 . 1 . . 

-Q^Oat, + — 0 ' 2 d a g + ^0:100,^036 


72 

2 2 
-CX-zdab + -C^lda/3 

2 . 1 . 

-0 3 0 a b + —Oild a g 

y io 


3 1 . 

TyO 2 + —OliOab 

lo 8 

3 1 . 

—a: 2 + -a:iOafe 
lo 8 


—a 1 
2 


Q, 

u, 

d, 
H, 
L, 

e. 


( 6 . 1 ) 


• fermion + anti-fermion —> gravitino + gaugino : The total cross-section for this 
channel is (l/32Mp) x (16a 3 + 7a 2 + llay/3), see row (I) of the table on page 
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[181] of the first reference in |15], For different degrees of freedom we then find: 


03 = 


—-— 6a X 

32M P 3 


&ab “I” 0^2 $a/3 T lrio ®l d a gS ab 


27 


108 

4 4 

-^Ol 3 5 ab + —«i5 a/ 3 

4 _ 1 _ 

-^jOi 3 d ab + —Oid a g 


1 1 r 

77^2 + T-zOtiOab 
O 

1 1 r 

-01-2 + — cMaft 

O 1Z 


-Q!l 

3 


< 3 , 

M, 

d , 

e. 


( 6 . 2 ) 


• sfermion + anti-sfermion —> gravitino + gaugino: The total cross-section for 
this channel is (1/32M|) x (8cc 3 + 7 ck 2 /2 + llai/6), see row (J) of the table 
on page [181] of the first reference in (Fj. For different degrees of freedom we 
then find: 


cr 4 = 


32 Ml 5ijX 


2 . 1 . 1 . . 
& 3 0 ab + — Ot20 a p + —— OLiO a pO ab 

lo 


27 


216 

2 2 
~^oi 3 5 ab + —CKif^g 

2 . 1 _ 
-^c>L 3 o ab + —aid Q/ 3 


+ - ai 6 ab 
1 1 r 

- 0.2 T - OiO nb 

16 24 


F" 1 

6 


< 3 , 

u, 

d, 

H, 

L, 

e. 


(6.3) 
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The total cross-section for the gravitino production will then be given by: 


Otot 


32M* SijX 


2 9 1 

+ —Oil dapfiab 


2 2 
-Oi^ab T yCqAa/3 


2 - 1 - 

0 CC3Cafe + QJiUq^ 
0 0 


9 3 \ 

— «2 + g^lhafe j 

9 3 . 

—do T — criOfih 
16 8 




< 5 , 

U, 

d, 

H, 

L, 

e. 


(6.4) 


As we pointed out in the previous section, particles and anti-particles associated to 
the bosonic and fermionic components of the multiplets which belong to an irre¬ 
ducible representation of a gauge group have the same parameter A r . This implies 
that 


3 2 3 

tot — ^ ^ ^ ^ ^ Ai t a,aAj t b,/3 (Otot^rel) 

i,j= 1 a,6=1 a, (3=1 


32M 2 ^ [6«3(2*4.q + A 2 U + A d ) + y a 2(3^.g + A 2 l + A 2 h ) 
+-oii(Aq + 8^ + 2A 2 d + 3.A| + 6A 2 e + 3 A 2 H )]. 


(6.5) 


The sum is taken over the three flavors of Q,u,d, L,e and the two Higgs doublets. 
After replacing E to t in Eq. (|A 6 |), and recalling that T max ~ m^/3, we obtain 



We remind that in full thermal equilibrium, E to t = ( 47 r/Mp) x (16ct3 + 6 ct 2 + 2a\) ~ 
( 10 _ 1 /Mp) (up to logarithmic corrections) B 0 . It is evident that the exact 
abundance of the gravitinos produced during the quasi-thermal phase depends on 
the composition of the reheat plasma. One expects the number of gravitinos thus 
produced to be maximum if the inflaton mainly decays to one flavor of LH (s)quarks 
(which are charged under the whole SM group). In this case Aq = 1/24 for the 


30 







relevant flavor 26 , while A = 0 for all other degrees of freedom. This results in a 
gravitino abundance: 


^ 3/2 


A 


3/4 / rri 

I -L max 

s v 228.75 ) V 10 10 GeV 


10 


-12 


(6.7) 


where A is given by Eq. (|5.1|) . 

An important point to note is that T max is accompanied by the factor A?^ in 
Eq. m- Therefore, despite the fact that T max « m^/3 can be as large as 10 12 GeV, 
the gravitino abundance can be at a safe level. First consider the case for unstable 
gravitinos. For T max ~ 10 12 GeV, the tightest bound from BBN (n 3 / 2 /s) < 1CV 16 
(arising for m 3 / 2 — 1 TeV and a hadronic branching ratio ~ 1) is satisfied if A < 10~ 6 . 
Much weaker bounds on A are found for a radiative decay. For example, A < 10 -3 
(1) if m 3/2 ~ 100 GeV (1 TeV). 

For stable gravitinos with C(keV) mass the dark matter limit, see Eq. 0) , is 
satisfied if A < 10 -6 . A much more relaxed bound A < 1 is obtained if m 3 / 2 — 100 
MeV. Therefore, in general, the gravitino production during the quasi-thermal phase 
is safe 27 . 

We conclude that late thermalization of the Universe due to SUSY flat directions 
eliminates the gravitino problem altogether in a natural way. 


7 . Leptogenesis in a quasi-thermal case 

7.1 Basic concept 

The baryon asymmetry of the Universe (BAU) parameterized as t/b — (^b — u b)/ s 
is determined to be 0.9 x 10 _1 ° by the recent analysis of WMAP data ||]. This 
number is also in good agreement with an independent determination from the pri¬ 
mordial abundance of light elements produced during BBN [Q. Any mechanism 
for generating a baryon asymmetry must satisfy: B— and/or L— violation, C— and 
CP—violation, and departure from thermal equilibrium 
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26 Tlie total number of degrees of freedom in one flavor of LH (s)quarks is 
2 (particle — antiparticle) x 2 (fermion — boson) x 2 (weak — isospin) x 3 (colour). 

2 'The above mentioned constraints on A are comfortably satisfied in a generic inflationary model 
where the inflaton is a gauge singlet and couples to the MS(SM) gravitationally, see the discussion 


in an Appendix 11.2, and in 11.3. 

28 Since B + L-violating sphaleron transitions are active at temperatures 100 GeV < T < 10 12 











Leptogenesis is an elegant mechanism which postulates the existence of RH neu¬ 
trinos, which are SM singlets, with a lepton number violating Majorana mass Mjy. It 
can be naturally embedded in models which explain the light neutrino masses via the 


see-saw mechanism J7DI . A lepton asymmetry can then be generated from the out-of¬ 
equilibrium decay of the RH neutrinos into Higgs bosons and light leptons, provided 
CP—violating phases exist in the neutrino Yukawa couplings [[71], [72], [73]]. The cre¬ 
ated lepton asymmetry will be converted into a baryonic asymmetry via sphaleron 
processes. 

In thermal leptogenesis the on-shell RH neutrinos whose decay is responsible for 
the lepton asymmetry are produced via their Yukawa interactions with the SM fields 
in a thermal bath |74j]. In SUSY there is RH sneutrinos which serve an additional 
source for leptogenesis [[75]]. This scenario works most comfortably if Tr > Mi > 
10 9 GeV @ § 29 . 

The decay of a RH (s)neutrino with mass Mj results in a lepton asymmetry 
via one-loop self-energy and vertex corrections [73]. If the asymmetry is mainly 


produced from the decay of the lightest RH states, and assuming hierarchical RH 
(s)neutrinos Mi M 2 ,M 3 , we will have ]57|] 

bn 3 — mi\ f Mi 


i|s~3x 10 


-10 


K 


0.05 eV ) V 10 9 GeV 


(full equilibrium), (7.1) 


for 0(1) CP-violating phases (mi < m 2 < m 3 are the masses of light mostly LH 
neutrinos). Here n is the efficiency factor accounting for the decay, inverse decay and 


scattering processes involving the RH states [38], [77 
A decay parameter K can be defined as 


K = 


Ti 

H(T = Mi) ’ 


(7.2) 


where Ti is the decay width of the lightest RH (s)neutrino. It can be related to an 
effective neutrino mass rhi such that K = mi/(ICC 3 eV), with the model-independent 
bound mi < rhi |5B| . 

GeV [6^], any mechanism for creating a baryon asymmetry at T > 100 GeV must create a B — L 
asymmetry. The final asymmetry is then given by B = a(B — L), where a = 28/79 in the case of 
SM and a = 8/23 for the MSSM [G9|] . 

29 There exist various scenarios of non-thermal leptogenesis |78], |?9|, |S(], |l[ |||, |S3|, |8l|, which 
can work for Tr < Mn- There are also leptogenesis models which implement soft SUSY breaking 
terms |[S6f . 
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If K < 1, corresponding to rhi < 1CT 3 eV, the decay of the RH states will be out 
of equilibrium at all times. In this case the RH states, which are mainly produced via 
scatterings of the LH (s)leptons off the top (s)quarks and electroweak gauge/gaugino 
fields |73|], never reach thermal equilibrium. The cross-section for producing the RH 
(s)neutrinos is oc T” 2 (M 2 ), when T > Mi (< Mi), and hence most of them are pro¬ 
duced when T ~ Mi. The efficiency factor reaches its maximum value for k ~ 0.1 
when rhi = 10 -3 eV. For larger values of rhi it drops again, because the inverse de¬ 
cays become important and suppress the generated asymmetry. Producing sufficient 
asymmetry then sets a lower bound, Mi > 10 9 GeV [38]. Successful thermal lepto- 
genesis therefore requires that Tr > 10 9 GeV. Note that this is at best marginally 
compatible with thermal gravitino production, see Eq. O)- 


7.2 Quasi-thermal leptogenesis 

The presence of flat directions slow down thermalization and lower the reheat tem¬ 
perature, such that, T r -C 10 9 GeV is naturally obtained (see Tables. 1 and 2). 
This is detrimental to thermal leptogenesis in the bulk of parameter space. Note 
that Eq. (O) implies that sufficient asymmetry will not be generated after the es¬ 
tablishment of a full equilibrium 30 . Needless to mention, like gravitino production, 
leptogenesis can still occur during the quasi-thermal phase and this is the topic of 
our interest in this subsection. 

In a quasi-thermal phase the reheat plasma is dilute, implying that the RH states 
are produced less abundantly than in full thermal equilibrium. Their abundance can 
be calculated from Eq. CT)> where the production cross-section is, a oc 1/M 2 . 
Note that the abundance depends on At, and hence on the composition of the reheat 
plasma. Significant production of the RH (s)neutrinos requires that the LH (s)leptons 
and/or the top (s)quarks be present in the reheat plasma. The relevant channels for 
producing the lightest RH neutrino, Ah, and sneutrino, Ah, are scatterings of LH 
(s)leptons with the largest Yukawa coupling to Ah, Ah off the LH top (s)quarks and 
anti-(s)quarks, and annihilation of top (s)quark-anti(s)quark pairs 31 . Summing over 

30 Thermal leptogenesis can work for Mi <C 10 9 GeV if the RH (s)neutrinos are degenerate |89| , 
or for specific neutrino mass models |j(|. 

31 The electroweak gauge/gaugino fields have a large mass and decay almost instantly, therefore, 
unlike the case with full equilibrium they do not participate in Ni , Ni production. 
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all processes (including weak-isospin and color indices) we find from Eq. (5.6) that: 


~ 10" 5 E 


Ah 


, n/r „ . /228.75\ 5/4 

(MpMi) (- 


A 


(7.3) 


where 


oc x 18 ( AlAq 3 + ALAt + AQ 3 At) , (7.4) 

is the production cross-section for Ni,Ni. Here Al, Aq 3 and A t denote the A 
parameter for the LH (s)leptons with the largest Yukawa coupling to N\, Ni, the LH 
top (s)quarks and the RH top (s)quarks respectively. Note that in full equilibrium 
they are all ~ 1. The three terms inside the parentheses are the contributions from 
the above mentioned processes respectively. The final baryon asymmetry generated 
in the quasi-thermal phase will then be given by: 

5/4 


Vb — 10 10 ( 228 /l ' 75 ) (A l Aq 3 + A L A t + Aq 3 A t ) k 1 ms mi 


A 


0.05 eV / V 10 9 GeV 


Mi 


(quasi — thermal) 


(7.5) 


If the inflaton mainly decays to the top (s)quarks, we have Aq 3 = At = A/ 36 32 , 
while A = 0 for all other degrees of freedom. This results in: 

3/4 


7/ b — 5 x 10 


-9 


K 


A 


Mi 


(7.6) 


228.75 J V 10® GeV. 

Since ?/b oc Mi, the maximum asymmetry is produced when Mi ~ 3 T max Ri3x 10 12 
GeV. For the largest efficiency factor k — 0.1, generating the correct asymmetry 
requires that A > 10~ 3 . This is compatible with the bound from the gravitino 
production for a radiatively decaying gravitino with m ^/2 — 100 GeV — 1 TeV (see 
the discussion in the previous section, Eq. 00- 

Note that from Eqs. 0) and ( fT5l) the gravitino abundance and the baryon 
asymmetry produced during a quasi-thermal phase are both oc A i A. However, due 
to the dependence on the composition of the reheat plasma, the marginality between 
leptogenesis and gravitino production can be very different from that in the case 
of full equilibrium. The best case scenario for leptogenesis occurs when the inflaton 
mainly decays to the top (s)quarks. In this case the marginality between leptogenesis, 
Eq. ( |7.6|) , and gravitino production, Eq. ( |6.7|) , is weakened by about one order of 
magnitude compared to the case of full equilibrium, see Eqs. O) and ( |7.1| ). 

32 Note that the total number of degrees of freedom in Q 3 and t is 36. 












There are other interesting differences which arise in the case of a quasi-thermal 
leptogenesis. Because the plasma is dilute in this phase, N\, Ni will not be brought 
into equilibrium even if ffi\ 1CU 3 eV. On the other hand, see Eq. the expan¬ 

sion rate of the Universe is (much) slower than the case with full thermal equilibrium 
when T ~ M\. This implies that out-of-equilibrium decay of N±, N± requires that 
rhi <C 10~ 3 eV. Moreover, since A L = 2 scatterings mediated by N\ , N\ are much 
less efficient in a dilute plasma, the resulting bound on Mi will be altered. 

To conclude, late thermalization of the Universe implies that thermal leptogenesis 
cannot generate sufficient asymmetry in the bulk of the parameter space. The new 
paradigm is the quasi-thermal leptogenesis, but this depends on the composition of 
the reheat plasma. If the inflaton mainly decays to the top (s)quarks, the marginality 
between gravitino production and leptogenesis will be improved (compared to the case 
in full equilibrium). As a result, the right amount of baryon asymmetry can be 
generated for sufficiently heavy RH (s)neutrinos. A more detailed and quantitative 
study is needed for a better understanding of the role of decays and inverse decays in 
a quasi-thermal leptogenesis. 


8. Thermalization after preheating 


So far we have focused on thermalization after the perturbative inflaton decay. Al¬ 
though for final reheating the perturbative decay is the most important, but we make 
some comments on the situation after a non-perturbative decay of the inflaton also. 
For detailed studies on thermalization after preheating see Refs. fU[ fl2[ [93|j . 


Consider a simple chaotic inflation model as in Appendix |11.1| with the following 
potential 

( 8 , 1 ) 


V ~ -mltf + fcW, 


2 j2, ,2 


where x is another scalar held. Here we have considered only the real parts of q 1 
and x- If h > 10 -6 , the inflaton oscillations decay to x quanta via broad parametric 
resonance 0. If h > 10 -4 , resonant production results in an extremely efficient 
transfer of energy from the zero-mode condensate in a typical time scale ~ lOOm^ 1 
(which depends weakly on the coupling h) g. 

Resonant particle production and re-scatterings lead to the formation of a plasma 
consisting of (j) and x quanta with typical energies 10 _1 (hm<pMp) 1/2 g. This 
plasma is in kinetic equilibrium but full thermal equilibrium is established over a 
much longer time scale than preheating i0 
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The occupation number of particles in the preheat plasma is 1 (which is op¬ 
posite to the situation after the perturbative decay). This implies that the number 
density of particles is larger than its value in full equilibrium, while the average 
energy of particles is smaller than the equilibrium value. It gives rise to large ef¬ 
fective masses for particles which, right after preheating, is similar to their typical 
momenta 0. Large occupation numbers also lead to important quantum effects due 
to identical particles and significant off-shell effects in the preheat plasma. Because 
of all these, a field theoretical study of thermalization is considerably more compli¬ 
cated in case of preheating. Due to the large occupation numbers, one can consider 
the problem as thermalization of classical fields at early stages [[Li], |92], [J8|| . In the 
course of evolution towards full equilibrium, however, the occupation numbers de¬ 
crease. Therefore a proper (non-equilibrium) quantum field theory treatment ]94 


will be inevitably required at late stages when occupation numbers are close to one. 

Similar complications also arise when considering particle production during 
thermalization. However let us make crude estimates based on Eqs. (|6.G| ) and (|7.5|) . 
At the end of preheating p 10 4 m^Mp and T ~ 10 1 (hrn^Mp)^ 2 . We therefore 
find from Eq. ( |5 . 1| ) that, A ~ h~ 2 228.75. Eq. ( |6.6| ) then implies overproduction 

of gravitinos in the preheat plasma unless, T max -C 10 10 GeV. This results in severe 
constraint on the models. For example, since T max > m^, it requires that m^ <C 10 10 
GeV which is a disaster from the point of view of inflaton generating the density 
perturbations. On the other hand, Eq. ( |7.5|) , implies that successful leptogenesis is 
now possible for M -C 10 9 GeV. In both the cases the situation is opposite to that 
after perturbative decay where del. 

We however caution the reader that these should be only taken as crude es¬ 
timates, since particle production from scatterings in the preheat plasma is more 
involved. In a densed plasma the scattering processes can be enhanced (for bosonic 
final states) or suppressed (for fermionic final states). Large effective masses can 
also kinematically suppress or shut-off some processes. We can nevertheless expect 
that our analysis captures the main qualitative aspects of particle production during 
thermalization after preheating. 


9. Conclusion 

In this paper we have presented a detailed account of the thermalization after infla¬ 
tion in SUSY and we discussed various implications. We have emphasized that the 
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final stage of reheating is the perturbative inflaton decay, even if the inflaton conden¬ 
sate decays non-perturbatively. For a wide range of inflaton couplings perturbative 
decay happens when the inflaton quanta dominate the energy density and therefore 
generates entropy. 

The most important result is the rate of thermalization is extremely slow in 
SUSY. It is often wrongly assumed that the inflaton decay products immediately 
thermalize if they have gauge interactions. Our message is that this (although true 
in realistic models in non-SUSY case) is not correct in SUSY. 

In any SUSY extension of the SM there is a large number of flat directions 
which are made up of squark, slcpton and Higgs fields. These flat directions acquire 
very large VEV during inflation which spontaneously break gauge symmetries in the 
early Universe. This induces very large masses to the gauge bosons (and gauginos) 
and suppresses main reactions which lead to kinetic and chemical equilibrium (i.e., 
the 2 —> 2 and 2 —> 3 scatterings with gauge boson exchange). As a result, the 
Universe enters a long period of a quasi-thermal phase during which the comoving 
number density and (average) energy of particles remain constant. This epoch lasts 
until the 2 —» 3 scatterings become efficient, at which point the number of particles 
increases and full equilibrium is established. The main results are given in Eqs. ( |2.13| ) 
and (fngp . 

Slow thermalization substantially lowers the reheat temperature of the Universe. 
The reheat temperature is practically decoupled from the inflaton decay and can be 
as low as O(TeV), even for large inflaton decay rates. It varies in a typical range, 
10 3 GeV < Tr < 10' GeV. This is underlined in Eq. (0 and demonstrated by the 
examples given in Tables. 1, 2. 


We studied particle production in a quasi-thermal phase. The general results are 
given in Eqs. (|5.2| ), (|5.3|) , (|5.6|) . Important aspect of a quasi-thermal particle pro¬ 
duction is its dependence on the composition of the reheat plasma (before complete 
thermalization). The abundance of particles thus produced does not depend solely 
on the maximum temperature in the quasi-thermal phase. We then specialized to 
two cases of physical interest, namely gravitino production and leptogenesis, with 
the results given in Eqs. and (|7.5|) respectively. 


The most important cosmological consequence of our study is the gravitino pro¬ 
duction. The Universe thermalizes at sufficiently low reheat temperatures which 
satisfy the tightest BBN bounds on thermal gravitino production. Our central mes- 
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sage is that there is a natural resolution to the infamous Gravitino problem lies 
within a consistent treatment of thermal history of the Universe within SUSY. We 
emphasize that the built in solution offered by SUSY renders any exotic modifications 
(such as late entropy release) unnecessary. Needless to say, this has very important 
implications for inflationary model building. 

On the other hand, quasi-thermal leptogenesis is necessary to generate sufficient 
baryon asymmetry since slow thermalization results in Tr <C 10 9 GeV, for which 
thermal leptogenesis does not work (unless in very special cases). As usual the ques¬ 
tion is the marginality between leptogenesis and gravitino production. Depending 
on the composition of the reheat plasma, this can be either relaxed or tightened 
compared to the case in full equilibrium. If the inflaton mainly decays into the top 
(s) quarks, it is possible to have successful quasi-thermal leptogenesis while keep¬ 
ing gravitino production under control. One can make further progress along these 
lines through more quantitative studies which takes into account of lepton-number 
violating interactions more carefully. 

To conclude, supersymmetry dramatically modifies thermal history of the Uni¬ 
verse, most importantly, it provides a built in solution in the form of flat directions 
which can naturally solve the gravitino problem. This, so far, neglected fact can 
remove one of the most serious obstacles for building consistent inflationary models 
in the framework of SUSY and in string inspired theory. 
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11. Appendix 


11.1 Last stage of inflaton decay 

Let us consider a simple model of chaotic inflation with a SUSY 33 superpotential 


W D + M?!'!', 


( 11 . 1 ) 


where $ is the inflaton superfield comprising of the inflaton 0 and the inflatino 0. 
It is coupled to another superfield T whose bosonic and fermionic components are 
denoted by y an d "0 respectively. Here we choose rn^ = 10 13 GeV, so that inflaton 
fluctuations generate the right amount of density perturbations. Eq. (|11.4| ) results 
in the following interaction terms in the scalar potential 


V D fcW + — =hm^<t>x 2 , 


( 11 . 2 ) 


where we have considered the real parts of 0 and x fields. A nice feature is that 
SUSY relates the couplings of the cubic 0y 2 and quartic 0 2 y 2 interaction terms of 
the inflaton. Note that the cubic term is required for complete decay of the inflaton 
field. 

At the end of inflation \<p\ ~ O(Mp). Preheating occurs if h > 10~ 6 , in which case 
the h 2 cj) 2 x 2 term takes over and, for h > 1CU 4 , leads to an explosive transfer of energy 
from the homogeneous condensate to x quanta ||. Eventually, after re-scattering 
of x quanta off the remaining condensate, a plasma is formed which consists of 
the same number 0 and y quanta with typical energies )$> which is in kinetic 
equilibrium |H], |92], |93|. This stage completes over a rather short period of time 


t ~ lOOm^ 1 ||. Full thermal equilibrium takes much longer to establish, but the 


33 A nice realization of chaotic inflation within supergravity is through implementing a shift sym¬ 
metry |)5). If inflaton Kahler potential has the form K = (0 + 0*) 2 /Mp, instead of the minimal 
form K = cf>*(f>/Mp, the scalar potential along the imaginary part of 0 remains flat even for Trans- 
planckian field values. Therefore it can play the role of inflaton in a chaotic model. Note that 
a shift symmetry also ensures that the (positive) Hubble induced corrections to the mass of flat 


directions vanishes at the tree-level as the cross terms in Eq. (4.2) disappear. It is also possible 
to realize chaotic inflation for sub Planckian field values in supergravity. For example, see the 
multi-axions fl9ql driven assisted inflation 1631. 













temperature of the resulting thermal bath will presumably be larger than my. This 
implies that the inflaton (and inflatino) quanta remain in thermal equilibrium as 
long as T > my 

Once T drops below my, due to Hubble expansion, the inflaton quanta become 
non-relativistic. The hm^cftx 2 /V 2 term then takes over, leading to a perturbative 
decay of the inflaton to (real and imaginary parts of) y, plus the fermionic partner 
0, at a rate Td = (h 2 / 8tt) m^. Note that regardless of how large h is, this stage of 
inflaton decay will be perturbative 34 . The reason is that, unlike the initial condensate 
there is no coherence among the decaying inflaton quanta at this stage. 

The inflaton quanta dominate the energy density of the Universe at the time of 
decay, and hence generate entropy, provided that 

9 

TTl 

n « jg. (ii.s) 


For my = 10 13 GeV, Eqs. (|1 1. 5|) and (|11.6| ) result in h < 1CU 2 . This is much 
weaker than the condition h < 10 -6 , which is required for the inflaton decay to 
be perturbative from the beginning. Therefore in SUSY a last stage of perturbative 
inflaton decay naturally follows preheating. Let us consider a simple model of chaotic 
inflation with a SUSY 35 superpotential 


W D -my$$ 




(11,4) 


where $ is the inflaton superfield comprising of the inflaton 0 and the inflatino 0. 
ft is coupled to another superfield T whose bosonic and fermionic components are 
denoted by y and 0 respectively. Here we choose my = 10 13 GeV, so that inflaton 
fluctuations generate the right amount of density perturbations. Eq. (|11.4|) results 


34 The situation is similar to that in the decay of SUSY partners of SM fields. These particles 
stay in thermal equilibrium at temperatures above their mass. Once T drops below their mass, 
they decay very quickly, but perturbatively through gauge couplings of 0(1). 

35 A nice realization of chaotic inflation within supergravity is through implementing a shift sym¬ 
metry |)5). If inflaton Kahler potential has the form K = (0 + 0*) 2 /Mp, instead of the minimal 
form K = ({)*(}>/Mp, the scalar potential along the imaginary part of 0 remains flat even for Trans- 
planckian field values. Therefore it can play the role of inflaton in a chaotic model. Note that 
a shift symmetry also ensures that the (positive) Hubble induced corrections to the mass of flat 
directions vanishes at the tree-level as the cross terms in Eq. (4.2) disappear. It is also possible 
to realize chaotic inflation for sub Planckian field values in supergravity. For example, see the 
multi-axions Ml driven assisted inflation l63|. 
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in the following interaction terms in the scalar potential 


V D h 2 (j) 2 x 2 + 


V2 


hm^x 2 , 


(11.5) 


where we have considered the real parts of (p and x fields. A nice feature is that 
SUSY relates the couplings of the cubic (px 2 and quartic <p 2 x 2 interaction terms of 
the inflaton. Note that the cubic term is required for complete decay of the inflaton 
field. 

At the end of inflation \(j)\ ~ O(Mp). Preheating occurs if h > 10~ 6 , in which case 
the h 2 (p 2 x 2 term takes over and, for h > 1 CU 4 , leads to an explosive transfer of energy 
from the homogeneous condensate to x quanta ||. Eventually, after re-scattering 
of x quanta off the remaining condensate, a plasma is formed which consists of 
the same number <p and x quanta with typical energies rn^ which is in kinetic 
equilibrium |)T], |92], |93|. This stage completes over a rather short period of time 
t ~ lOOm ^ 1 [|J. Full thermal equilibrium takes much longer to establish, but the 
temperature of the resulting thermal bath will presumably be larger than m^. This 
implies that the inflaton (and inflatino) quanta remain in thermal equilibrium as 
long as T > m 0 . 

Once T drops below due to Hubble expansion, the inflaton quanta become 
non-relativistic. The hm^cpx 2 /\/2 term then takes over, leading to a perturbative 
decay of the inflaton to (real and imaginary parts of) x, plus the fermionic partner 
ip, at a rate Td = (h 2 / 87 t) m<y Note that regardless of how large h is, this stage of 
inflaton decay will be perturbative 36 . The reason is that, unlike the initial condensate 
there is no coherence among the decaying inflaton quanta at this stage. 

The inflaton quanta dominate the energy density of the Universe at the time of 
decay, and hence generate entropy, provided that 


T ‘«W 


( 11 . 6 ) 


For m<£ = 10 13 GeV, Eqs. (|11.5|) and (|11.6| ) result in h < 1CT 2 . This is much 
weaker than the condition h < 10 -6 , which is required for the inflaton decay to 
be perturbative from the beginning. Therefore in SUSY a last stage of perturbative 
inflaton decay naturally follows preheating. 


36 The situation is similar to that in the decay of SUSY partners of SM fields. These particles 
stay in thermal equilibrium at temperatures above their mass. Once T drops below their mass, 
they decay very quickly, but perturbatively through gauge couplings of 0(1). 











11.2 Gravitationally decaying inflaton 


37 As a first example, we consider a model of Ref. |97]] in which the inflaton sector is 
gravitationally coupled to the MSSM sector. The scalar potential in supergravity is 


given by f4l[ 


V = e a l GiG' - ] Mp, 


where G is the Kahler function and in minimal supergravity is defined as 


G = it + log 


\wf 

Ml 


(11.7) 


( 11 , 8 ) 


The Xi denote the scalar helds in the theory, and lower and upper indices on G 
denote its derivative with respect to Xa and x* respectively. The inflaton sector 
superpotential and the MSSM superpotential Wmssm are given by 


W* = (d> - Mp ) 2 


Umssm = Uijk^i^ j^k- 


(11.9) 


Here $ and T,; denote the inflaton and the MSSM chiral superfields, respectively, and 
Uijk are the MSSM Yukawa couplings. The minimum of inflaton potential is located 
at 0 = Mp, around which it takes the form m 2 (0 — Mp) 2 . This model realizes new 
inflation in minimal supergravity. Obtaining density perturbations of the correct 
size from quantum fluctuations of the inflaton requires that m<p = 10 13 GeV |)7 
Eq. (|11,71) leads to the following term 


VijkJT^XiXjXk, ( 11 . 10 ) 

1V1 p 

in the scalar potential. If m$ is much larger than the soft SUSY breaking scalar 
masses, the partial width for inflaton decay to three scalars is ~ y 2 m^/M|. For 
m<p = 10 13 GeV this is always the case, particularly in models with weak scale 
SUSY. 

There is also a term 


e G/2 [G ij - G i G j - G k G i l] 0^-, 


( 11 . 11 ) 


in the Lagrangian [[U]] which describes the couplings of fermionic partners of T, 
denoted by x, in the two-component notation 


yijk-g-Axii’jA ■ 


( 11 . 12 ) 


37 This subsection is motivated by our discussion with Antonio Masiero. 
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It results in a partial width for inflaton decay to two fermions and one scalar which 
is same as that for decay to three scalars ~ y 2 m^/Mp. This implies that inflaton 
decay produces the same number of particles and sparticlcs. Because of the large 
top Yukawa coupling y t « 1, the inflaton in this model mainly decays to the top 
(s)quarks, LH bottom (s)quarks, Higgs H u and Higgsino H u . The total inflaton 
decay rate is therefore 

T d ~ 10 ~ 2 ml/M 2 . (11.13) 


Since the inflaton decays into three-body final states, and hence the average energy 
of decay products is (E) ~ m$/ 3, and hence T max ~ m^/9. Eq. (|5TD then implies 
that 


A 


m^Y 

) 


(11.14) 


For = 10 13 GeV, Eq. ( |11.14| ) leads to A ~ 10 -11 . This implies that the reheat 
plasma is extremely dilute, and hence substantially far from full equilibrium. For 
different degrees of freedom we have: 


Ai 




(11.15) 


where yi denotes the superpotential Yukawa coupling of the i —th degree of freedom. 
Hence it is the largest for top (s)quarks, LH bottom (s)quarks, H u and H u . 


11.3 Right-handed sneutrino as the inflaton 

It is also possible that the inflaton is directly coupled to some of the MSSM fields. 
This happens, for example, in the model of Refs. |)8], |99| where one of the the RH 
sneutrinos N plays the role of the inflaton. The relevant part of the superpotential 
in this case reads: 

W D \m n NN + hiHJALi , (11.16) 

where N is the multiplet containing the RH sneutrino N which plays the role of the 
inflaton (and its fermionic partner N), and hi are the Yukawa couplings governing 
the inflaton decay. With an appropriate choice of non-minimal Kahler function, the 
scalar potential remains flat at large held values TV > Mp, and this model realizes 
chaotic inflation in supergravity |98|, (see also ||100|| ). Quantum fluctuations of the 
sneutrino result in density perturbations of the correct size, provided that M N — 10 13 
GeV. 
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The inflaton in this case mainly decays into the LH (s)leptons, H u and H u . Note 
that the same number of particles and sparticles are produced in inflaton decay, so 
long as Mat is much larger than soft SUSY breaking masses (which is the case for 
weak scale SUSY). The total inflaton decay rate is then given by 


r d 


if 

47T 


M 


N 



(11.17) 


Since the inflaton decays into two-body final states, we have (E) = m^/2 right after 
the decay completes, implying that T max ~ m^/Q. Eq. ( |5.1| ) then results in 

For Mm = 10 13 GeV and 10~ 6 < h < 10” 3 we hnd 10“ 12 < A < 1. For the i —th 
(s) lepton doublet we have 


The (s)lepton singlets, (s)quarks, gauge fields and gauginos are not produced in two- 
body decays of the inflaton. However they are inevitably produced at higher orders 
of perturbation theory |5(J, therefore, they have much smaller but non-vanishing 
values of A. 


11.4 Flat directions and thermalization: additional considerations 

The main reason behind slow thermalization of the Universe in SUSY is that the flat 
direction VEV, and hence the mass of gauge bosons, remains large for a sufficiently 
long time. In Section |2] we have assumed that the flat direction VEV is just redshifted 
by the Hubble expansion while oscillating. Here we consider further details of the 
flat direction dynamics and examine their importance. 

• Flat direction decay: 

In our analysis we have assumed that flat directions do not decay until the 
Universe fully thermalizes. The flat directions have gauge and Yukawa cou¬ 
plings, generally denoted by y , to other fields which results in a decay rate 
= (y 2 /47t) m 0 . Note, however, that the flat direction VEV induces a mass 
y\tp\ for the decay products. The decay is therefore kinematically forbidden 
until y\(p\ < m 0 . On the other hand, the flat direction decay to particles to 
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which it is not directly coupled is kinematically allowed at all times. However, 
such decays are mediated by the fields which are coupled to the flat direction, 
and hence suppressed by a factor ~ ( mo/y\(p\ ) 2 relative to the leading order 
decay. 


In both cases it turns out that the flat direction decays at a time when the 
expansion rate is given by: 


< 


1 rrin 


A , . ( 11 . 20 ) 

4?r M 

The flat direction will not decay before the establishment of full thermal equi¬ 
librium if H ^ < i? t hr where, depending on the case, H t h r is given by Eqs. ( j2.13| ) 
and (12.15|) . This is generically the case, in particular for the examples of Ta¬ 
ble. 1. 


One might also wonder whether the flat direction could promptly decay via 
preheating since y\(p\ m o at the onset of its oscillations ||. However, SUSY 
breaking soft mass terms in general lead to out-of-phase oscillations of the real 


and imaginary parts of the flat direction [45]. As shown in [101], a tiny effect 
of this type is sufficient to shut-off resonant decay of the flat direction. Similar 
effects will be present if two or more flat directions with non-zero VEV oscillate 
slightly out of phase. 


• Scatterings off the flat direction condensate: 

Energetic particles in the reheat plasma scatter off the zero-mode quanta in 
the flat direction condensate. In order for the flat direction to affect ther- 
malization, it should survive against evaporation by such scatterings (at least) 
until the Universe fully thermalizes. Scatterings which are mediated by gauge 
and gaugino fields, similar to diagrams in Fig. (1), play the main role here. 
However, the center-of-mass energy in such scatterings is s 1//2 ~ (4 Ethq) 1 ^ 2 . It 
turns out that for H > H t h r we typically have AEm 0 -C a]^] 2 . The evaporation 
rate is therefore given by 


r eV a ~ 


10a 


nEmo 

~W' 


( 11 . 21 ) 


which is < F t h r . This implies that the flat direction evaporates only after the 
establishment of full equilibrium. 


Note that the condensate contains zero-mode quanta with a number density 
n,„ = m 0 \ip\ 2 . For sufficiently large values of (fo that affect thermalization, 
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is much larger than the number density of particles in the reheat plasma n. 
These zero-mode quanta can participate in the 2^3 scattering similar to 
those shown in Fig. (2). Then one might wonder whether scatterings off the 
flat direction condensate would result in a thermalization rate larger than what 
we obtained in Eq. (|2.10|) . 

However emitting a gauge boson with a mass ~ g\<p\ requires that 4 Em 0 > 
a|<p| 2 (note that {AEmf) 1 ^ 2 is the center-of-mass energy). As just mentioned, 
we have the opposite inequality for H > H t h r - This implies that the 2 —> 3 
scatterings off the condensate do not have enough energy to produce on-shell 
gauge bosons. Inelastic scatterings which happen at higher orders (mediated by 
off-shell gauge and gaugino fields) can increase the number of particles but their 
cross-section will be suppressed by (at least) a factor of 4Em 0 /a|(^| 2 relative 
to those in Fig. (2). This results in a rate 

r thr - 10« 2 ^, (11-22) 

m 


where we have used n v = m 0 \ip\ 2 . This turns out to be subdominant to r t hr, 
see Eq. (|2.10|) , whenever the flat direction VEV is large enough to affect ther¬ 
malization. 


• Early oscillations of the flat direction: 

The flat direction starts oscillating once the Hubble expansion rate drops below 
its mass, which we have considered to be the soft breaking mass m q. 


The thermal effects can trigger early oscillations of the flat directions ||102| . |103| . 
This follows from a general consideration that the fields which are coupled 
to a flat direction (once kinematically accessible to a thermal bath) induce 
a plasma mass for it. If the plasma mass exceeds the expansion rate of the 
Universe at early times, it will trigger early oscillations of the flat direction. 
When the reheat plasma is in full thermal equilibrium, this indeed happens 
in many cases [|102| , |103|| . Earlier oscillations of the flat direction also imply 
earlier thermalization of the reheat plasma, see Eqs. (|2.13| ) and (|2.15| ), and 
hence higher reheat temperatures. 


Before the establishment of full equilibrium, however, the reheat plasma is 
dilute and the resulting plasma mass ~ (n/E) 1 ^ 2 is smaller by a factor of 
(Ul/228.75) 1 ^ 4 . Moreover, particles which are coupled to the flat direction have 
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a large mass, unless these couplings are very small, and hence decay via gauge 
interactions almost instantly. Note that the inverse decays are inefficient in a 
dilute plasma. Therefore the plasma masses are completely negligible and will 
not affect the flat direction dynamics. 
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